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Abstract 

We study the Q 2 variation of the first moment of the nucleon's spin-dependent 
structure function G\. As Q 2 — > the moment is determined by the low en- 
ergy theorem for Compton scattering. In the deep-inelastic region the moment 
is calculated using twist expansion to order 1/Q 2 - Based on these limits, we 
construct a formula which smoothly interpolates between the two regions. 
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Recently, polarized deep-inelastic scattering has proven to be an excellent tool for study- 
ing the spin structure physics of the nucleon |]-|| . Supplemented with the operator product 
expansion analysis in Quantum Chromodynamics (QCD), experimental data at high-energy 
provides a direct measurement of the matrix elements of spin-dependent operators in the 
nucleon. A much discussed example in the current literature is the axial charge, or the for- 
ward matrix element of axial current, whose measurement by the EMC collaboration casts 
doubt on our traditional understanding of the proton's spin structure M. 

A closely related question is can one learn anything about the nucleon's spin structure 
from electro-production experiments away from the deep-inelastic limit? In particular, what 
insight do the spin structure functions G\ and G2 provide at low and moderate Q 2 1 Not 
long ago, Anselmino et al. Q pointed out that at the real photon point (Q 2 = 0), the first 
moment of G\ (called the sum rule in the following text) is related, via the celebrated Drell- 
Hearn-Gerasimov (DHG) sum rule ||, to the anomalous magnetic moment of the nucleon, 
and thus the physics of the G\ structure function again appears simple in the Q 2 — > limit. 
Together with knowledge from the deep-inelastic limit, the authors in Ref. j| constructed a 
model for the sum rule at all Q 2 . This has motivated a number of proposals to measure G\ 
and G2 at low energy ||,[7|. 

In Ref. ||, one of us pointed out that the analysis made in Ref. H excluded the nucleon's 
elastic contribution to the moment, which in the Q 2 —>■ limit dominates the entire inelastic 
contribution calculated from the DHG sum rule. He argued that the moment has to include 
this contribution if it is to be analyzed in twist expansion in the deep-inelastic limit and 
its experimental measurement is to be used to extract the matrix elements of higher-twist 
operators. 

In this Letter we study the Q 2 variation of the sum rule by exploring the physics of the 
Q 2 — > and Q 2 — > 00 limits. In the first limit, we rely on the low energy result derived in 
Ref. ||, to calculate the exact value and the first derivative of the sum rule at Q 2 = 0. In 
the second limit, we use a twist expansion appropriate for the deep inelastic region, focusing 
on the 1/Q 2 correction term. The matrix elements of higher twist operators are related to 
moments of the quark distributions functions gif(x)(i = 1,T, 3 and / = u,d, s...) through 
a novel use of the QCD equations of motion, which are in turn evaluated in the MIT bag 
model. As an application, we discuss the correction of higher twists to the Bjorken sum rule 
in the deep inelastic limit. Having obtained analytical results valid for the low and high 
ranges of Q 2 , we construct a simple parameterization to smoothly interpolate both limits, 
which should be checked experimentally. 

To begin we consider the following fixed-mass sum rule, 

Q 2 f°° ^ / „ 2 \^ 



where Gi(u, Q 2 ) is one of the nucleon's spin dependent structure functions in the nucleon 
tensor, 

W* = -ie^ p q a [S^ + ^(uS* - P?(S ■ q))]. (2) 

The lower integration limit in Eq. (JJ) implies the elastic contribution to G\ is also included. 
Here P and S are the nucleon's momentum and polarization, q is the virtual photon mo- 
mentum, M is the nucleon mass, v = P ■ q and Q 2 = —q 2 (e 0123 = 1). In deep-inelastic 
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limit, one defines scaling functions gi(x,Q 2 ) = vjM 2 G\ and g2{x,Q 2 ) = (u/M 2 ) 2 G2- The 
sum rule then becomes, 

T(Q 2 )= C g x {x,Q 2 )dx, (3) 
Jo 

which is just the first moment of the scaling function. 

Let us first consider the small Q 2 behavior of the sum rule. Introduce a spin-dependent 
virtual-photon Compton amplitude Si(u, Q 2 ) whose imaginary part is proportional to G\, 
and write down the unsubtracted dispersion relation, 

5i(^ 2 ) = 4f -#^GiKg 2 ). (4) 

JQ 2 /2 V ,Jl — V 1 

Through this, we relate the sum rule to the Compton amplitude at v — 0, 

r(Q 2 ) = ^i(o,Q 2 ). (5) 

At small v and Q 2 , the dominant contribution to S\ comes from the nucleon pole diagrams 

a 

Q>) = + F 2 )[^±- - ^ - Fl (6) 

where F\ and F 2 are the Dirac and Pauli form factors of the nucleon. From this, we obtain 
for Q 2 -> 0, 

nQ 2 ) = \F 1 (F 1 + F 2 )~^ f2 F 2 Q 2 . (7) 

This result can be shown to be accurate up to the order of Q 2 in the small Q 2 region by 
explicitly evaluating Eq. ([!]): The elastic contribution to Gi is proportional to 5{2u — Q 2 ), 
producing the first term in Eq. ([?]); the integral over inelastic contributions is just the DHG 
sum rule in the limit of Q 2 — > and the second term in Eq. (|7|) reproduces this in the same 
limit. 

The elastic contribution vanishes identically at Q 2 = because of energy-momentum con- 
servation, and the DHG sum rule indicates T(0) = 0. Thus, due to the elastic contribution, 
T(Q 2 ) is non-analytic around Q 2 = 0, i.e., 

r(Q 2 = o) + r(Q 2 -> o). (8) 

To remedy this, one can take two approaches: The first approach subtracts away the elastic 
contribution from the sum rule for Q 2 ^ 0. The new sum, T(Q 2 ) = T(Q 2 ) — 1/2F 1 (F 1 + F2)), 
is a smooth extension of the DHG sum rule to virtual-photon scattering. The approach we 
take in this paper is to redefine r at Q 2 = 0, 

T(Q 2 = 0)=T(Q 2 ^0). (9) 

This approach ensures that the sum rule at low Q 2 can be treated with the twist expansion 
that we will discuss below. The expansion is for moments of the gi structure function which 
include the integration limit x = 1, where the elastic contribution resides. 
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Since Eq. (|7|) is accurate up to the order of Q 2 , we can determine T(Q 2 ) and its first 
derivative in Q 2 — > limit, 





r p (o) 




r n (o) 


dTP(Q 2 


), 


dQ 2 


-|q 2 =o 


dT n (Q 2 


), 


dQ 2 


-\Q 2 =0 



-8.631GeV" 2 , 

-0.479GeV- 2 , (10) 

where p, n refer to proton or neutron and the squares of the proton and neutron charge 
radii (r 2 ) c . r . = (0.862 fm) 2 and (r 2 ) c . r . = —(0.342 fm) 2 have been used. The initial slope 
of T P (Q 2 ) is primarily determined by the elastic contribution as the inelastic contribution, 
-Kp/8M 2 = — 0.455GeV~ 2 , is only about 5% of the total. Therefore, one expects that for 
small Q 2 , T P (Q 2 ) is mainly given by the elastic contribution. In contrast, due to a numerical 
coincidence, the elastic part of T n (Q 2 ) is negligible compared with the inelastic part. 

In the limit of large Q 2 (Q 2 >>> Aq CD ), T(Q 2 ) can be calculated in terms of the twist 
expansion, 

r=2,4... [Q ) 2 

where /i r (Q 2 ) are matrix elements of quark-gluon operators which scale like Aqq 2 -,. The in- 
dependence in [l t are logarithmic and can be calculated in perturbative QCD. If the nucleon 
mass were zero, [i T (Q 2 ) would contain only twist-r operators. The effect of the nucleon mass 
is to induce contributions to fi T {Q 2 ) from lower twist operators, as we shall illustrate below. 
The leading term in Eq. flTT| ) is well-known, 

A*2 = 2 E e / fl o/' ( 12 ) 

f=u,d,s... 

where the summation covers quarks of all flavors / and ao/ is the axial charge defined by 
the matrix element of axial current = ^fj^js^f. (PS\A^\PS) = 2aofS tl . The QCD 
radiative corrections have been calculated to the first order in a s (Q 2 ) for the singlet 
contribution (af = 2(aQ U + aod + ao s )/9) and to the third order ]10| for the non-singlet 
contribution (a^ 3 = (2a 0u — aod — a os)/9). The proton-neutron difference of the moment 
defines the Bjorken sum rule, 



1-^2^) -3.58 1 I -20.2 
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(13) 



where g& = a 0u — a 0c z = 1.257 is the neutron decay constant. 

The 1/Q 2 power corrections to T were first studied by Shuryak and Vainshtein (SV) 
IIJ . Using the collinear expansion technique ]T5[, one of us has calculated in Ref. |13] the 



entire 1/Q 2 corrections to the g\ scaling function in terms of a few multi-parton distribution 
functions. Specializing to the first moment, we find 
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^ = q E 4 h/ + rf 2/ - 4/ 2/ ] m 2 
y / 

= (A + L> + F)M 2 . (14) 

where A = \Yjf e } a 2f comes from the twist-two contribution and a 2 / is defined as, 

(P3\i>fy^j s iD' tl iD^ f \P3) = 2a 2f S {u P^P^\ (15) 

with (■ ■ •) denotes symmetrizing the indices and subtracting the trace; D — | J^/ e /^2/ comes 
from the twist-three contribution and d 2 f is defined as, 

(PSlgTpF^^iplPS) = 2d 2f S [a P^ l] P^\ (16) 

with [• • •] denotes anti-symmetrizing the indices and F afl1 = l/2e atJlial3 F a p is the dual of 
the gluon field tensor; F = — |Z)/ e //2/ comes from the twist-four contribution and f 2 f is 
defined as, 

(PS\gi> f F^^ f \PS) = 2f 2f M 2 S». (17) 

We note that the result quoted for fi 4 in Ref. |IJ is (2 A + 2D + F)M 2 . 

To study the QCD radiative corrections to /14, one has to consider operator-mixing from 
gluon operators, the anomalous dimensions of which are not currently available and their 
matrix elements are difficult to estimate. Therefore, in the following discussion, we neglect 
entirely the scale dependence of /Z4. 

The higher twist operators in Eqs. (|16|) and ([[7]) depend explicitly on gauge fields. To 
calculate their matrix elements we need a wave function of the nucleon containing gluon 
components. However, for special types of higher twist operators such as the present case, 
we can eliminate the gluons in terms of the "bad" components of quark fields using the QCD 
equations of motion [|13j. Then the higher twist matrix elements can be related to moments 
of parton distributions with no explicit gluon fields. Indeed, by defining in the light-cone 
gauge (A • n = 0), 

<7(i,T,3)/(*) = \j d ^(PS\4, f Q W) 7 ^ f (\n)\PS), (18) 

where Q\ = n, Qt = —St/M, Q 3 = —2p/M 2 and n and p are two null vectors (n 2 = p 2 = 
and p ■ n — 1), we find, 

1 r ,, 



d 2 f = - J x (3g 2f (x) + 2g lf (x))dx, 

hf = \j x 2 (7gif(x) + 1202/0*0 - 9g 3f (x))dx, (19) 

where g 2 f — 9rf — 9if- These relations are exact in QCD. 

We choose to estimate the 1/Q 2 corrections to the sum rule in the simplest version 

of the MIT bag model, in which the bag boundary simulates gluon confinement ]T4] , |r5 . 

Using gif(x) (i = 1,T, 3) calculated in this model, we obtain for the proton, A p = 0.0065, 
Dp = 0.0092, and F p = 0.0155. Inserting them into Eq. <M), we have, 
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/4(Bag) = 0.031M 2 . 



(20) 



Compared with the size of \F 2 = 0.126 ± 0.025 from the EMC data or 0.175 from the Ellis- 
Jaffe sum rule, the bag l/Q 2 power correction is about 10% at Q 2 = 2 GeV 2 and about 
2% at Q 2 = 10 GeV 2 . Assuming there are no abnormal twist-six or higher contributions, 
we conclude that most Q 2 variations of the proton sum rule occur below 1 GeV 2 . For the 
neutron, the bag model predicts, A n = D n = F n = 0, and the l/Q 2 correction vanishes: 

/^(Bag) = 0. (21) 

This follows from the SU(6) structure of the bag wave function, which also predicts fi 2 = 0. 

The higher twist matrix elements have also been calculated using the QCD sum rule 
(QSR) technique by Balitsky et al ||16|| . Their most recent result in terms of our notation is 

/4(QSR) = -(0.023 ± 0.015)M 2 , 

/^(QSR) = -(0.006 ± 0.004)M 2 . (22) 

Thus, the QCD sum rule calculation gives a power correction the same size as the bag 
calculation with the correction for the neutron being significantly smaller than for the proton. 
However, the sign of the correction differs from the bag result in Eq. (|20|) . This difference 
has a large effect on the Bjorken sum rule at small Q 2 . 

The Bjorken sum rule has recently been extracted from the data on the proton |l| and 
neutron [pf gi structure functions, 

f 1 gf~ n (x, 2 GeV 2 )dx = 0.146 ± 0.021 [2], 
Jo 

= 0.152 ±0.025 [T?j. (23) 

In QCD, the Bjorken sum rule at low Q 2 is contaminated by higher twist corrections dis- 
cussed above. If the QSR result is used for the correction, one obtains a theoretical prediction 
at the same Q 2 , 

f 1 gl~ n {x, 2 GeV 2 )dx = 0.160 [17], (24) 
Jo 

On the other hand, the bag result produces, 

f 1 g\~ n {x, 2 GeV 2 )dx = 0.182. (25) 
Jo 

While Eq. ( |2"4"D gives a corrected Bjorken sum rule within the experiment errors, the bag 
calculation disagrees with the extraction of the sum rule by the E142 collaboration by 1.7cr 
and with the extraction by Ellis-Karliner by 1.2a. In our opinion, a deviation from the 
Bjorken sum rule means either a measurement of higher twist matrix elements, or that the 
data are inconsistent. 

From the high and low Q 2 knowledge of the sum rule, we propose a model for T P (Q 2 ) in 
the entire Q 2 region, 

!-((,») = 1^^0(1-^)^11^, (26, 
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where the first term is the elastic contribution with its derivative modified by the Q 2 term. 
The second term is basically a+b/ Q 2 and the denominator serves to suppress the contribution 
at small Q 2 . From the EMC data and the various constraints derived above, we determine 
all Aj except A 4 , which controls the size of the twist-six contribution. The solid and upper- 
dashed curves shown in Fig. 1 are our parameterization with the bag and QSR higher twist 
matrix elements, respectively. [We choose A4 = 0.3, which gives a n§ ~ —0.03.] The dotted 
curve represents the result of the twist expansion to order 1/Q 2 and the dot-dashed curve 
represents the elastic contribution. As can be seen from the figure, the different choices for 
higher twist matrix elements result in about 15% difference in T in the Q 2 = 0.5 to 1.0 
GeV 2 region. A similar interpolation is made for the neutron, and the result is shown as the 
lower-dashed curve. 

Thus it appears that the Q 2 variation of the T(Q 2 ) sum rule is quite simple. Nevertheless, 
its experimental measurement is interesting, particularly around Q 2 = 0.5 GeV 2 . If we know 
T(Q 2 ) in an extended Q 2 region, we can fit data with a parameterization similar to the one 
used in Eq. (|26|). Then by expanding in a 1/Q 2 power series, we can extract the higher- twist 
matrix elements, such as f2f, which shall provide valuable insight into the spin structure of 
the nucleon. 

We thank Bob Jaffe for checking Eq. (IT3) and I. Balitsky and A. Vainshtein for several 
useful discussions. PU acknowledges the support of a NSERC scholarship from the Canadian 
government. 
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FIGURE CAPTION 

Fig. 1: A model for the sum rule T(Q 2 ) at all Q 2 . 
The meaning of the curves is explained in the text. 
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